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Abstract 

The Lamb shift in muonic hydrogen continues to be a subject of experimental and 
theoretical investigation. Here my older work on the subject is updated to provide a 
complementary calculation of the energies of the 2p-2s transitions in muonic hydrogen. 

Introduction 

The energy levels of muonic atoms are very sensitive to effects of quantum electrodynamics 
(QED), nuclear structure, and recoil, since the muon is about 206 times heavier than the 
electron pQ. In view of a proposed measurement of the Lamb shift im muonic hydrogen [2], 
an improved theoretical analysis seems to be desirable. Since the first theoretical analysis 
[3j, the subject of the Lamb shift (the 2p-2s transition) in light muonic atoms has been 
investigated with increasing precision by a number of authors [H El 13 El 13 UD] • The 
present paper provides an independent recalculation of some of the most important effects, 
including hyperfine structure, and a new calculation of some terms that were omitted 
in the most recent literature, such as the virtual Delbriick effect fTP. An alternative 
calculation of the relativistic recoil correction is presented. 

In the numerical calculations the fundamental constants from the CODATA 1998 ( 12 ) 
are used, i.e.: He, m^, rUe, m„ = 137.0359998, 197.32696 MeV-fm, 105.658357 MeV, 

0.5109989 MeV, 931.4940 MeV, respectively. The changes in these constants in the CO- 
DATA 2002 compared with CODATA 1998 are too small to make any relevant difference 
in the results. 

Vacuum Polarization 

The most important QED effect for muonic atoms is the virtual production and anni- 
hilation of a single e'^e~ pair It has as a consequence an effective interaction of order 
aZa which is usually called the Uehling potential (^3J This interaction describes 
the most important modification of Coulomb's law. Numerically it is so important that it 
should not be treated using perturbation theory; instead the Uehling potential should be 
added to the nuclear electrostatic potential before solving the Dirac equation. However, 
a perturbative treatment is also useful in the case of very light atoms, such as hydrogen. 
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Unlike some other authors, we prefer to use relativistic (Dirac) wave functions to 
describe the muonic orbit. This is more exact, and as will be seen below, it makes a 
difference for at least the most important contributions. The wave functions are given in 
the book of Akhiezer and Berestetskii ^3] and will not be given here. In perturbation 
theory, the energy shift due to an effective potential AV is given by 

1 /"OO /"OO 

AE„« = — ■/ q'dqAViq)- drMqr)[Fl + GU (1) 

Jo Jo 



where Fnn and Guk are the small and large components of the wave function, n is the 

i and +i if j = i - i. 



principle quantum number and k is equal to — + 1) if j = i + h and +i if j — ^ ^ 



AV{q) is the Fourier transform of the physical potential. 
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AV{q) = An- I ■ jo(gr) • AV{r) dr (2) 
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AV^(r) = ^ - / q'-Uqr)-AViq)dq (3) 







As is well-known the Uehling potential in momentum space is given by 

Vuehiiq) = -i^Ml . GE{q) ■ F{<P) = -4n{aZ) ■ GE^q) ■ 
where Ge is the proton charge form factor, sinh((/)) = q/{2me) and 

= ]- + (coth2(0) - 3) ■ [1 + ■ coth(0)] (4) 

U2{q) is defined in The vacuum polarization corrections were calculated in momentum 
space; the formulas (124,125,127) of ^ are completely equivalent to (200) in [TU]. If 
the correction to the transition 2pi/2 — 2si/2 is calculated in lowest order perturbation 
theory using nonrelativistic point Coulomb wave functions, the result is 205.0074 meV, in 
agreement with other authors ^U] . 

The same procedure was used to calculate the two-loop corrections; the corresponding 
diagrams were first calculated by Kallen and Sabry |Tnj. The Fourier transform of the 
corresponding potential is given in 011]. The result for a point nucleus is 1.5080 meV. 

In momentum space including the effect of nuclear size on the Uehling potential is 
trivial, since the corresponding expression for AV{q) is simply multiplied by the form 
factor. The numbers obtained were the same for a dipole form factor and for a Gaussian 
form factor, provided the parameters were adjusted to reproduce the experimental rms 
radius of the proton. The correction can be regarded as taking into account the effect 
of finite nuclear size on the virtual electron-positron pair in the loop. The contribution 
of the Uehling potential to the 2p-2s transition is reduced by 0.0081 meV with a proton 
radius of 0.862 fm ^7j, and by 0.0085 meV with a proton radius of 0.880 fm ^Hj. This 
result is consistent with the number given in pOJ (eq.(266)). More recent values for the 
proton radius have been given by Sick [T^ (0.895 ± 0.018 fm) and in the newest CODATA 
compilation I2DI (0.875 ±0.007fm). 
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The numerical values given below were calculated as the expectation value of the 
Uehling potential using point-Coulomb Dirac wave functions with reduced mass: 





point nucleus 


i?p=0.875fm 


Uehling 
Kaellen-Sabry 


2j9i/2 - 2Si/2 2p3/2 - 2Si/2 

205.0282 205.0332 
1.50814 1.50818 


2pi/2 - 2Si/2 2p3/2 - 2Si/2 

205.0199 205.0250 
1.50807 1.50811 



The effect of finite proton size calculated here can be parametrized as -0.0109(r^). 
However higher iterations can change these results. For a very crude estimate, one can 
scale previous results for helium |3] and assume that the ratio of nonperturbative to 
perturbative contributions was the same, giving a contribution of 0.175meV. 
The contribution due to two and three iterations have been calculated by |H] and |2S1, 
respectively, giving a total of 0.151meV. An additional higher iteration including finite 
size and vacuum polarization is given in ref. |8j (equations(66) and (67)) and ref. [TU] 
(equations(264) and (268)). These amount to -0.0164(r^). The best way to calculate this 
would be an accurate numerical solution of the Dirac equation in the combined Coulomb- 
plus Uehhng potential. 

The mixed muon-electron vacuum polarization correction was recalculated and gave 
the same result as obtained previously, namely 0.00007 meV. PT| ITU]. 

The Wichmann-KroU |22] contribution was calculated using the parametrization for 
the potential given in [T]. The result obtained (-0.00103 meV) is consistent with that 
given in jjljj, but not with that given in jS]. 

The equivalent potential for the virtual Delbriick effect was recomputed from the 
Fourier transform given in JT] and pp. The resulting potential was checked by reproducing 
previously calculated results for the 2s-2p transition in muonic helium, and the 3d-2p 
transitions in muonic Mg and Si. The result for hydrogen is -|-(0. 00135 ± 0.00015) meV. 
As in the case of muonic helium, this contribution very nearly cancels the Wichmann-KroU 
contribution. The contribution corresponding to three photons to the muon and one to 
the proton should be analogous to the light by light contribution to the muon anomalous 
moment; to my knowledge, the corresponding contribution to the muon form factor has 
never been calculated. It will be comparable to the other light by light contributions. For 
an estimate, the correction to the Lamb shift due to the contribution to the anomalous 
magnetic moment was calculated; it amounts to (-)0. 00002 meV; the contribution to the 
muon form factor is one of the most significant unknown corrections. 

The sixth order vacuum polarization corrections to the Lamb shift in muonic hydrogen 
have been calculated by Kinoshita and Nio j23j. Their result for the 2p-2s transition is 

AE^^) = 0.120045- (aZ)2 -m,. (^-j^ ^ 0.00761 meV 

It is entirely possible that the as-yet uncalculated light by light contribution will give a 
comparable contribution. 

The hadronic vacuum polarization contribution has been estimated by a number of 
authors [211 123 CHI- It amounts to about 0.012 meV. One point that should not be 
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forgotten about the hadronic VP correction is the fact that the sum rule or dispersion 
relation that everyone (including myself) used does not take into account the fact that 
the proton (nucleus) can in principle interact strongly with the hadrons in the virtual 
hadron loop. This is irrelevant for the anomalous magnetic moment but probably not for 
muonic atoms. An estimation of this effect appears to be extremely difficult, and could 
easily change the correction by up to 50%. Eides et al. [TU] point out that the graph 
related to hadronic vacuum polarization can also contriibute to the measured value of the 
nuclear charge distribution (and polarizability) . It is not easy to determine where the 
contribution should be assigned. 



Finite nuclear size and nuclear polarization 

The main contribution due to finite nuclear size has been given analytically to order {aZy 
by Friar [211. The main result is 



2aZ ( aZrrir 



n 



I 1\ OiZTXlrf , O, 



(2) 



+ {aZf{FREL + mlFMR) 



(5) 



where (r^) is the mean square radius of the proton. For muonic hydrogen, the coef- 
ficient of (r^) is 5.1975 (meV fm~^), giving an energy shift (for the leading term) of 
3.862±0.108 meV if the proton rms radius is 0.862±0.012) fm. The shift is 4.163±0.188 meV 
if the proton rms radius is 0.895±0.018) fm, and 3.979±0.076 meV if the proton rms radius 
of 0.875±0.007) fm. . The second term in Eq. (0) contributes -0.0232 meV for a dipole form 
factor and -0.0212 meV for a Gaussian form factor. The parameters were fitted to the pro- 
ton rms radius. This can be written as -0.0347(r2)3/2 or 0.0317(r2)3/2^ respectively. This 
differs slightly from the value given by Pachucki [0]. The model dependence introduces 
an uncertainty of about ±0.002 meV. The remaining terms contribute 0.00046 meV. This 
estimate includes all of the terms given in [211, while other authors [Hj give only some 
of them. Clearly the neglected terms are not negligible. There is also a contribution of 
-3 • 10~®meV to the binding energy of the 2pi/2-level, and a recoil correction of 0.012 meV 
to the binding energy of the 2s-level. 

As mentioned previously, the finite-size contributions to vacuum polarization can be 
parametrized as - 0.0109(r2) - 0.0164(r2), giving a total of -0.0273(r2) or -0.0209(6) meV 
if the proton radius is 0.875 fm. 

The contribution due to nuclear polarization has been calculated by Rosenfelder [2Z] 
to be 0.017 ± 0.004 meV, and by Pachuki ^ to be 0.012 ± 0.002 meV. Other calculations 
[2H1 I2n] give intermediate values (0.013 meV and 0.016 meV, respectively). The value 
appearing in table [T] is an average of the three most recent values, with the largest quoted 
uncertainty, which is probably underestimated. 



Relativistic Recoil 

As is well-known, the center-of-mass motion can be separated exactly from the relative 
motion only in the nonrelativistic limit. Relativistic corrections have been studied by 
many authors, and will not be reviewed here. The relativistic recoil corrections summa- 
rized in include the effect of finite nuclear size to leading order in m^/mjv properly. 
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Up to now this method has been used to treat recoil corrections to vacuum polarization 
only in the context of extensive numerical calculations that include the Uehling potential 
in the complete potential, as described in p. They can be included explicitly, as a 
perturbation correction to point-Coulomb values. Recall that (to leading order in l/m^v), 
the energy levels are given by 

E = E,-^ + -^{hir) + 2BoPiir)) (6) 

where Er is the energy level calculated using the reduced mass and Bq is the unperturbed 
binding energy. Also 

h{r) = -P,{r)[P,{r) + -Q2{r)] - ^Q,{r)[P,{r) + Q,{r)/r'] (7) 

r Sr 

Here 

POO 

Pi{r) = AiraZ r' p{r')dr' = -V{r) - rV'{r) (8) 

J r 

Q2{r) = AnaZ I r'^p{r')dr' =rV(r) 
Jo 

Jo 

An effective charge density pvp for vacuum polarization can be derived from the 
Fourier transform of the Uehling potential. Recall that (for a point nucleus) 

aZ 2a 

Vuehi[r) = — ^— •Xi(2mer) 

where Xn{x) is defined in |T]. In momentum space, the Fourier transform of 'W'^V is 
obtained by multiplying the Fourier transform of V by —q^. Note that using the nor- 
malizations of P^Ej, one has V^V" = — 47raZp where p is the charge density. One then 
obtains 



2a p (^2 _ 1)1/2 ^ 1 ^ ^2 r <l'-U<ir) , 



2 

TT 
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(9) 



q U2{q)jo{qr) dq 



U2{q) is defined in [T]. It is also easy to show that 
dVuehi aZ2a 



dr r 371 



-Xi(2mer) + 2meXo(2mer) 
r 



1 2a 2me 

— Vuehi{r) + [aZ)- Xo(2mer) 

r Stt r 



Keeping only the Coulomb and Ueliling potentials, one finds 



Pi{r) 
Q2{r) 



-aZ^(2me)xo(2mer) 

OTT 

2a 

aZ { 1 + — [xi(2mer) + (2mer)xo(2mer)] 



aZ 
2 

71 



2a 
3^ 



dz- 



1)1/2 



1 

2? 



[6qr — {qr)^] cos(gr) + [3(gr)^ — 6] sin(gr) 



+ Am'lz'^ 



dq 



where Xn{x) is defined in jTj. Corrections due to finite nuclear size can be included when 
a model for the charge distribution is given. This done by Friar [211 (and confirmed 
independently for two different model charge distributions); the contribution due to finite 
nuclear size to the recoil correction for the binding energy of the 2s-level is -0.013 meV. The 
factor l/m„ is replaced by l/(m^ + ttin), also consistent with the calculations presented 
in 

Since vacuum polarization is assumed to be a relatively small correction to the Coulomb 
potential, it will be sufficient to approximate Q2{'r) by aZ/r. After some algebra, one can 
reduce the expectation values to single integrals: 



2a 

■ 2mi,aZ — 
37r 



{z' - 1)1/2 



az 



az 



(1 



azj 



1 + az)^ 



On 



(10) 



dz 



r 



I 2a; 
[aZ) m,.me — 

OTT 



2{az\ 



{Z^ - 1)1/2 



2(1 + azY 



5m + 



1 



2(1 + azY 



1 

' 2^2 
5ii 1 dz 
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with a = 2me/{aZmr). When Eq. (jlOj) is multiplied by —2BQ/{m^ + m^) this results 
in a shift of — 0.00015 meV for the 2s-state and of — 0.00001 meV for the 2p-state, and 
when Eq. (II Ij) is multiplied by l/{m^ + mj^) this results in a shift of 0.00489 meV for the 
2s-state and of 0.00017 meV for the 2p-state. These expectation values also appear when 
vacuum polarization is included in the Breit equation jST]. 
Finally, 



aZ , 



aZYml2a r {z^-lfl^ 



37r 



6 

az 
3 + 

az) 



1 

2^2 



2 + a^ 2 , 
ln(l + az^ 

1 + az az 



3(azV + 2az-l 



(1 



az) 



+ 



(12) 



4(1 + 



1 - 3a;z - 2{az)^ , 
"^^0 H TTTT- rA ^^1 



4(1 + 



az, 



dz 
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When multiplied by l/{m^ + mi^) this results in a shift of 0.002475 meV for the 2s-state 
and of 0.000238 meV for the 2p-state. 

Combining these expectation values according to equations IH] and [71 one finds a 
contribution to the 2p-2s transition of -0.00419 meV. To obtain the full relativistic and 
recoil corrections, one must add the difference between the expectation values of the 
Uehling potential calculated with relativistic and nonrelativistic wave functions, giving 
a total correction of 0.0166 meV. This is in fairly good agreement with the correction 
of .0169 meV calculated by Veitia and Pachucki [3T], using a generalization of the Breit 
equation [321 which is similar to that given in [U]. The treatment presented here has the 
advantage of avoiding second order perturbation theory. 

The review by Eides et.al gives a better version of the two photon recoil (Eq. 136) 
than was available for the review by Borie and G. Rinker pp. Evaluating this expression 
for muonic hydrogen gives a contribution of -0.04497 meV to the 2p-2s transition. Higher 
order radiative recoil corrections give an additional contribution of -0.0096 meV p!"0| . 
However, some of the contributions to the expressions given in |TU] involve logarithms 
of the mass ratio m^/mN- Logarithms can only arise in integrations in the region from 

to m-Ar; in this region the effect of the nuclear form factor should not be neglected. 
Pachucki [S] has estimated a finite size correction to this of about 0.02 meV, which seems 
to be similar to the term proportional to {r^)(2) given in Eq.© as calculated in the 
external field approximation by Friar This two-photon correction requires further 

investigation. In particular, the parametrization of the form factors used in any calculation 
should reproduce the correct proton radius. 

An additional recoil correction for states with £ 7^ has been given by (see also 
[ID]). It is 

---^(-«(;^) '«) 

When evaluated for the 2p-states of muonic hydrogen, one finds a contribution to the 
2p-2s transition energy of 0.0575 meV for the 2pi/2 state and -0.0287meV for the 2p3/2 
state. 



Muon Lamb Shift 

For the calculation of muon self-energy and vacuum polarization, the lowest order (one- 
loop approximation) contribution is well-known, at least in perturbation theory. Including 
also muon vacuum polarization (0.0168 me V) and an extra term of order {ZaY as given 
in [lU]: 




which contributes -0.00443 meV, one finds a contribution of -0.66788 meV for the 2si/2 — 2pi/2 
transition and -0.65031 meV for the 2si/2 — 2^3/2 transition. 

A misprint in the evaluation of the contribution of the higher order muon form factors 
(contributing to the fourth order terms) has been corrected. The extra electron loop 
contribution to -F2(0) is should be 1.09426(a/7r)^. This reproduces the correct coefficient 
of (a/vr)^ from the muon (g-2) analyses. This is .7658, which is equal to 1.09426-0.32848. 
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The fourth order electron loops [30] dominate the fourth order contribution (-0.00169 meV 
and -0.00164 meV, respectively). The rest is the same as for the electron pp. The 
contribution of the electron loops alone is -0.00168 meV for the 2si/2 — 2pi/2 transition 
and -0.00159 meV for the 2si/2 — 2^3/2 transition. 

Pachuki has estimated an additional contribution of -0.005 meV for a contribution 
corresponding to a vacuum polarization insert in the external photon. 



Summary of contributions 

Using the fundamental constants from the CODATA 1998 (^21) one finds the transition 
energies in meV in table ^ Here the main vacuum polarization contributions are given 
for a point nucleus, using the Dirac equation with reduced mass. Some uncertainties have 
been increased from the values given by the authors, as discussed in the text. 

The finite size corrections up to order (aZ)^ can be parametrized as 
5.1975(r2) - 0.0109(r2) - 0.0164(r2) + 0.0347(r3)(2). 

In the case of the muon Lamb shift, the numbers in tabled are for the 2si/2 — 2pi/2 
transition. The corresponding numbers for the 2si/2 — 2^3/2 transition are -0.65031 meV 
and -0.00164 meV, respectively. 

Fine structure of the 2p state 

There are two possible ways to calculate the fine structure. One is to start with the point 
Dirac value, include the contribution due to vacuum polarization, as calculated above, as 
well as the spin-orbit splitting (computed perturbatively) due to the muon's anomalous 
magnetic moment, and recoil as given by Eq. p3|) . The results are summarized in table El 

An alternative method is to use the formalism given in [Hj (and elsewhere, see, eg. 
jS21lin|) which gives the energy shift as the expectation value of 

I dV 1 + + (a^ + l/2)mAr/m^ ^ 

— ■ L ■ cr^ [14:) 

r ar yriNm^ 

Note that 

1111 



rrii^imfj, 2m^ 2m^ 2m^ 

so that the terms not involving in the spin-orbit contribution are really the Dirac fine 
structure plus the Barker-Glover correction (Eq. IT^ 

The Uehling potential has to be included in the potential V{r). For states with i > 
in light atoms, and neglecting the effect of finite nuclear size, we may take 



IdV _ aZ 
r dr 



2a jz' - ly/' f 1 
37r J, V 



1 + ^ / ^ ■ ( 1 + ) • (1 + 2merz) ■ dz 



(15) 



which is obtained from the Uehling potential |T31 by differentiation. Then, assuming 
that it is sufficient to use nonrelativistic point Coulomb wave functions for the 2p state, 
one finds 

1 \ / 1 \ 
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Contribution 


Value (meV) 


Uncertainty (meV) 


Uehling 


oner riooo 




Kallen-babry 


l.oUol 




W icnmann-Kroil 


-0.00103 




virt. Delbrueck 


0.00135 


0.00015 


mixed mu-e Vr 


0.00007 




hadronic VP 


U.Uii 


U.UUz 


sixtli order 1231 


0.00761 




recoil PO] (eql36) 


-0.04497 




recoil, higher order |lUj 


-U.UUyO 




recoli, nnite size j/bj 


0.013 


r\ r\r\ i 

0.001 


recoil correction to Vr [Ij 


-0.0041 




additional recoil [32j 


0.0575 




muon Lamb shiit 






second order 


-0.66788 




fourth order 


-U.UUiby 




nuclear size (i?p=0.875fm) 




0.007 fm 


main correction |2H| 


-3.979 


0.076 


order {aZf ^ 
order [ctZ)^ ^26] 


0.0232 


0.002 


-U.UUUo 




correction to VP 


-0.0083 




polarization 


0.015 


0.004 


Other (not checked) 

VP iterations jH] 

VP insertion in self energy jS] 

additional size for VP [lOj 


0.151 
-0.005 
-0.0128 





Table 1: Contributions to the muonic hydrogen Lamb shift. The proton radius is taken 
from (20) . 





- E{2pi,2) (meV) 


Dirac 


8.41564 


Uehling(VP) 


0.0050 


Kallen-Sabry 


0.00004 


anomalous moment 




second order 


0.01757 


higher orders 


0.00007 


Recoil (Eq.jni) 


-0.0862 


Total Fine Structure 


8.352 



Table 2: Contributions to the fine structure of the 2p-state in muonic hydrogen. 
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where 

2a {z^ - 1)1/2 / 1 \ / 1 2az \ 

~ 3^ A i2 y^2^)' \{l + azf ^ {l + azf J 

with a = 2mf./{aZmr). The result for the fine structure is 

-{aZfml _ / 1 

n^{2£+l)K \miymfj_ 2m^ m^m, 

where is given by Eq. lfTBj) . In this case, the terms involving in the expression for 
the muon Lamb shift are included, and should not be double counted. With a numerical 
value of = 0.000365, one finds a contribution of 0.00305 meV (compared with 0.005 meV 
using Dirac wave functions). 

Numerically, the terms not involving give a contribution of 8.3291 meV and the 
contribution from gives a contribution of 0.0176 meV, for a total of 8.3467 meV, in 
good agreement with Eq. 80 of [8^ . When the vacuum polarization correction is added, the 
result is only very slightly different from the Dirac value of 8.352 meV. The contribution 
due to the anomalous magnetic moment of the muon is the same in both cases. 

In both cases one should include the i?2/2Mjv-type correction to the fine structure, 
(see jnij Eq(38)). This is tiny (5.7 ■ 10~^meV) and is not included in the table. Friar 
has given expressions for the energy shifts of the 2p-states due to finite nuclear size. 
These were calculated and found to give a negligible contribution (3.1 • 10~^meV) to the 
fine structure of the 2p-state. 




Hyperfine structure 

The hyperfine structure is calculated in the same way as was done in earlier work [HI EI, 
but with improved accuracy. Most of the formalism and results are similar to those given 
by 0. 



The 2p state: 

The hyperfine structure of the 2p-state is given by [B] {F is the total angular momentum 
of the state) 



/ldV\ ,^ 
- — ) ■ l + K 



Am^m AT \ r dr / 2p 



2(1 + x)5,y(F(F+l)- 11/4) 



+ 6jy(C,.i(l + a^)-2(l + a;)) , , 



F 1 



2 2 



1 F 1 

2 2"' 



:i8) 



where j = a/2j~+T, the 6-j symbols are defined in [SS] , and 
Cpi = 5fi - 2Sfo - (1/5)(5f2- 



X 



mf,{l + 2Kp) 



2mjv(l + Kp) 

represents a recoil correction due to Thomas precession The same correction due 

to vacuum polarization (Eq. (fTHjl ) should be applied to the HFS shifts of the 2p-states, as 
well as to the spin-orbit term. 
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As has been known for a long time [3 13 E] , the states with total angular momentum 
F = 1 are a superposition of the states with j = 1/2 and j = 3/2. Let the fine structure 
sphtting be denoted hj 6 = £'2^3/2 — -E'2pi/2, and let 

and 13' = (3- (1 + e2p)- 

The energy shifts of the 2p-states with total angular momentum F (notation "^^^^Lj) 
are then given in table El 

State Energy Energy in meV 

'^]hi'2 -p'{2 + x + a^)/8 -5.971 

^j9i/2 (A - R)/2 1.846 

3p3/2 (A + R)/2 6.376 

^P3/2 5 + /3'(l + 5x/4-a^/4)/20 9.624 

Table 3: Hyperfine structure of the 2p-state in muonic hydrogen. 



where 

A = 6-(3'{x-a^)/16 
= [S- (3'{l + 7x/8 + a^/8)/6]2 + {f3')\l + 2x- a^f/289> 

(Here 5 = 8.352 me V) Some minor errors in have been corrected. These numbers differ 
slightly from those given in ref. fTT?] . 

The 2s-state: 

The basic hyperfine splitting of the 2s-state is given by 

Aup = ^^^^ ■ (1 + Kp) ■ (1 + du) =/?■(! + «m) = 22.8332 mel/ 
3m^m7v 

(see, for example [TU] , Eq. (271,277)) As was shown in [Hllin], the energy shift of the 
2s-state is given by: 

AEas = /5 ■ (1 + a^) ■ (1 + Syp + e^^rtex + ^Breit + ^FS.rec) " \^F\ " 35fo]/4 (19) 

Here (jSU) 

l''(a^)' , ,„-4 



''Vertex 



and (W, Eq. (277)) 

^Breit = 7; = 1-13 ■ 10 
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The vacuum polarization correction has two contributions. One of these is a result of a 
modification of the magnetic interaction between the muon and the nucleus and is given 
by (see jTj) 



£VPl 



4« r . fR„Jr] 
r ar 



37r2 



/?n.(0) 



1 + 



Q^jo{qr)GA.i{q) dq 



dz 



(20) 



2z'^ J Am'^[z'^ + {q/2me 



One can do two of the integrals analytically and obtains for the 2s-state (with a = 2me/ {aZrUr 
and sinh(0) = g/(2me) = K/a) 



^VPl 



4a 
3^ 



;F{(P)GM{aZmrK) dK 



7 



+ 



6 



(21) 



where F{(f)) is known from the Fourier transform of the Uehling potential and is given by 
EqQ. 

The other contribution, as discussed by EH] arises from the fact that the lower 
energy hyperfine state, being more tightly bound, has a higher probability of being in a 
region where vacuum polarization is large. This results in an additional energy shift of 



Vuehi {r)iJ2s {r)6Mi^2s {r)d^r 



Following Ref. [SI] with y = {aZmr/2) ■ r, one has 



aZnir 



exp(-y) 



(l-y)(ln(2y)+7) + 



131/ - 3 - 2i/2 



Ay 



(7 is Euler's constant), and 

'^2s{r) = V2^(0)(l - 2/)exp(-|/) 
One finds after a lengthy integration 

evP2 = ^ r T^^GE{aZmrK)F{<P) 



37r2 7o l + K'^ 



17 



41 



24 



+ 



2 (1 + 7^2)2 {i + K^f (1 + 7^2)4 
ln(l + ir2; 



1 + 7^2 



7 



+ 



6 



+ 



K 



1 - 



19 20 

+ 



12 



2(1 + 7^2) (1 + 7^2)2 + 



(22) 



SternheimjSS] denotes the two contributions by 6m and 6e, respectively. An alternative 
exression, obtained by assuming a point nucleus, using Eq.(131) from [1 for the Uehling 
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potential, and doing the integrations in a different order, is 



£VP2 — 




371 J ^ ■ V J ■ (1 + azy 



az 1 23 3 



(23) 



2 1 + az 8(l + az)2 2(1 + a^)^ 



ln(l + az) ■ (l h — — - — 777 1 



dz 



with a = 2me/ {aZmred)- Both methods give the same result. 

In the case of ordinary hydrogen, each of these contributes "ia^ = 1.997 • 10^^. The 
accuracy of the numerical integration was checked by reproducing these results. One can 
thus expect that muonic vacuum polarization will contribute 3a'^/A ^ 4 ■ 10~^, as in the 
case of normal hydrogen. This amounts to an energy shift of 0.0009 meV. Contributions 
due to the weak interaction or hadronic vacuum polarization should be even smaller. 
For muonic hydrogen, one obtains £ypi=0. 00211 and £:yp2=0. 00325 for a point nucleus. 
Including the effect of the proton size (with Ge{(i) = Gm{(i) as a dipole form factor) 
reduces these numbers to 0.00206 and 0.00321, respectively. For the case of muonic ^He 
[7j, the corresponding numbers are eypi=0. 00286 and £yP2=0. 00476. The contribution 
to the hyperfine splitting of the 2s-state is then 0.0470 meV+0. 0733 me V=0. 1203 meV 
(0.1212 meV if muonic vacuum polarization is included). The combined Breit and vertex 
corrections reduce this value to 0.1207meV. (0.1226meV if the proton form factors are 
not taken into account). 

The contribution to the hyperfine structure from the two loop diagrams JH] can be 
calculated by replacing U2{aZmrK) = (a/37r)F(0) by U4{aZmrK) (as given in [Bll]) in 
equations l^andl^ The resulting contributions are 1.64-10"^ and 2.46-10"^, respectively, 
giving a total shift of 0.0009 meV. 

The correction due to finite size and recoil have been given in as -0.145 meV, while 
a value of -0.152 meV is given in jHHJ. Ref. jH] also gives a correction as calculated by 
Zemach (|33) equal to -0.183 me V, but claims that this correction does not treat recoil 
properly. The Zemach correction is equal to 



where {r)(2) is given in [Hll^niEZl- Using the value (r)(2) = 1.086 ± 0.012fm from |37i| , 
gives ezem = —0.00702, and a contribution of of -0.1742 meV to the hyperfine splitting 
of the 2s state. Including this, but not other recoil corrections to the hyperfine structure 
of the 2s-state gives a total splitting of 22.7806 meV. Additional higher order corrections 
calculated in Ref. ^38j amount to a total of -0.0003 meV and are not included here. 

Summary of contributions and Conclusions 

The most important contributions to the Lamb shift in muonic hydrogen, including 
hyperfine structure, have been independently recalculated. A new calculation of some 
terms that were omitted in the most recent literature, such as the virtual Delbriick effect 
[TT] and an alternative calculation of the relativistic recoil correction have been presented. 



^ Zem 



2aZmr{r){2) 
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Numerically the results given in table ^ add up to a total correction of 
(206.032(6) - 5.225 (r2) + 0.0347 (r2)3/2) nieV = 202.055±0. 12 meV. (for the value of the 
proton radius from [ID])- As is well known, most of the uncertainty arises from the 
uncertainty in the proton radius. 

However, the contribution of the light-by-light graph to the muon form factor has not 
yet been calculated. Also, since m^/mp = 0.1126 is much larger than aZ, it is possible 
that recoil corrections of higher order in the mass ratio, that have never been calculated, 
could be significant at the level of the expected experimental accuracy of about 0.01 meV. 
In particular, the two-photon recoil corrections, including finite nuclear size, should be 
recalculated to resolve (small) inconsistencies among various theoretical results. 
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